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I. Summary 

A class of explicit and implicit total variation diminishing (TVD) schemes for the compressible 
Filler and Navier-Stokes equations has been developed 11-4]. They have the property of not 
generating spurious oscillations across shocks and contact discontinuities. In general, shocks 
can be captured within 1-2 grid points. For the inviscid case, one can divide these schemes into 
upwind TVD schemes and symmetric (non-upwind) TVD schemes. The upwind TVD scheme is 
based on the second-order TVD scheme developed by llarten (5). The symmetric TVD scheme 
developed by the author is a generalization of Roe [6] and Davis’s [7] TVD Lax-WendrofT scheme. 
The objective of this paper is to investigate the performance of those schemes on some viscous 
and inviscid airfoil steady-state calculations. A comparison of the symmetric and upwind TVD 
schemes is included. 


II. Description of Algorithm for System of Hyperbolic Conservation Laws 

The notion of upwind and symmetric TVD schemes, including formulation and extension to 
system cases (in uniform Cartesian grids), can be found in references (1-3,5). Here the extension 
of the implicit second-order-accurate TVD scheme for hyperbolic systems of conservation laws 
in curvilinear coordinates (4) is briefly described. 

Consider a two-dimensional system of hyperbolic conservation laws 


dQ dl'(Q) 0G(Q) = 

Ot dx dy 


( 2 . 1 ) 


Here Q , F(Q) and G(Q) are column vectors of m components. 

A generalized coordinate transformation of the form £ -• £(x,y) and t; — »;(x,y) which 
maintains the strong conservation law form of equation (2.1) is given by 


w + *Ml + ?C(5J = 0> ( 2 . 2 ) 

at at a,i 

where Q ~ Q/J> F -■ {ZxF I £ V G)/J i G — (r} x F + ij v G)/J, and J -- Car r /v ~ Uie Jacobian 
transformation. Let A -■ dF/dQ and B = 8G/dQ\ then the Jacobians A and B of F and G 
can be written as 

A=tf*A-K„2>) ( 2 - 3a ) 

B - (ii z A + r} v B). (2.3b) 

Let the eigenvalues of A be (nl,o|, ...,a7*) and the eigenvalues of B be (o* ,a*, ...,a”). Denote 
and R n as the matrices whose columns arc eigenvectors of A and and denote R^ and 
/f" 1 as the inverses of R$ and R n . 

Let the grid spacing be denoted by A£ and Ary such that £ ~ jAf and tj — A:Ar/. Let Q } ^ i k 
denote some symmetric average of Q 3t k and Q 3 + i t k (for example, ~ 0.5 + + i,fc + Qji*)’ 
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or the Roc’s average (8) for gas dynamics). Let + ^, ll } + j, ll j *, denote the quantity.- a‘ { , li(, 
7^ 1 related to A evaluated at Q yl i k . Similarly, let a' fc + ,, /i fc + a. ^ k 'i denote the quantities 
a[,, /i„, /f„ l related to /i evaluated at 4?, |fc+ j. 


Define 


1 j , k_ 

J + * i+i6.5'*'("j J4 T,* W J(k j 


(2.4a) 


as the dilference of the characteristic variables in the local (•direction, and define 


o,., i -- R 


... i _ 

fc+ 3 6.5 * (j],k+l -f J),k) 


(2.4b) 


as the difference of the characteristic variables in the local rj-direction. The J ]>k is the Jacobian 
transformation evaluated at. (j'Af.fcArj). The averaged Jacobians are used here in order to 
preserve frccstrcain. 

With the above notation, a one-parameter family of TVD schemes can be written as 

h a'* < a-* 

- < 5 ?.* A«n - *>[f;, - a*(i - ♦) •- g;.». ,] • (**») 

where 0 < 0 < 1, A* - At/ A£, = At/Ari , and At is the time step. This onc-paramcter 

family of schemes contains implicit as well as explicit schemes. When 0 ~ 0, (2.5) is an explicit 
method; when 0 / 0, it is an implicit scheme. For example, if 0 • = 1/2, the time differencing 
is the trapezoidal formula and scheme (2.5) is second-order in time and space. If 0 “ 1, the 
time differencing is the backward Euler method and scheme (2.5) is first-order in time but 
second-order in space. 

The numerical flux function i k for both the upwind and symmetric. £VD schemes [2,3] 
can be expressed as 

Fyi i k = - d* Fj-t i,it •+■ * (2.5b) 

Similarly, we can define the numerical flux G j k + 1 . 

Upwind TVD Scheme [5] 

For a particular form of the upwind TVD scheme [2,5], the elements of the $ J+ i denoted by 
(<A; + i)-,f = are 


(^‘ + i) y - ^(ai + |)(4 + ^+i) + '» < i4iK-+* 

3' - ininmod(a;._j ,a‘ + i). 


with 


(2.Ca) 


(2.6b) 


The minmod function of a list of arguments is equal to the smallest number in absolute value if 
the list of arguments is of the same sign, or is equal to zero if any argument is of opposite sign. 
The function V ; * s 

*w={[*J«+.*)/!. |:| <: • < 2<ic > 


Here 1 is a small positive parameter (see reference 5 for a formula of e), and 


d+i ” (0 o* + [ = 0, 


(2.6d) 
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where a 1 . 4 are elements of (2.4a). 

j t 8 

Symmetric TVD Scheme |3| 

For two particular forms of the symmetric TVD scheme, the elements of the 4 > 3+ ji denoted 

l>y Wj,art 1 

(2,a) 

and , 

j )* S ---• 1 4 ) [a'lH " minmod(a'^,tt' + ., 1 a^ i )] . (2.71.) 

A Conservative Linearized ADI Form For Steady-State Application 

For two-dimensional steady-state applications, the implicit schemes (2.5) can be solved by 
some appropriate relaxation method other than alternating direction implicit (ADI) and will 
be the direction of future research. The schemes considered in this paper are implemented in 
a conservative noniterative ADI form | 2 |. For steady-state applications, the numerical solution 
is independent of the time step. The implicit operator has a regular block tridiagonal structure 
and the resulting block tridiagonal matrix is diagonally dominant. One can modify a standard 
central difference code by simply changing the conventional numerical dissipation term into the 
one designed for the TVD scheme; i.e., the third term of equation (2.5b). -The only difference 
in computation is that the current scheme requires a more elaborate dissipation term for the 
explicit operator; no extra computation is required for the implicit operator. For the Navicr- 
Stokes applications, the convection terms are discretized by a TVD scheme, and the diffusion 
terms are discretized by central difference approximations. 

A conservative linearized ADI form of equation (2.5) used mainly for steady-st.at.e applications 
as described in details in reference ( 2 ) can be written as 


' ' ***«/. 1 .. “ 


D' ~ -A 4 


F n , 

J + Q I* 


Vi.* 


A" 


■in 


>.*4 k 


<V..i 


where 


= Q n q D 

~ 2 /* + ».*+ n 5+*A. 


// 




n 


(3.1a) 

(3.1b) 

(3.1c) 

(3.Id) 
(3.1e) 


The nonstandard notation 




(3. If) 


is used and 


n >+i.fc D ‘ = diag 

n l k U fl = dia 8 


- max rj>{a l }+ ,) 


- max V>(aj [+ ,.) 


uw-*>;.*> 

(D Jlk+1 - D,, k ). 


(3‘g) 

(3.1h) 


Here A j4Iik , H Jik+l are ( 2 . 3 ) evaluated at (j + 1 ,k) and {j,k + 1 ). The expression diag(i‘) 
denotes a diagonal matrix with diagonal elements z l . All of the inviscid calculations shown in 
this paper use (3.1). _ 
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For steady-state application, a simple algorithm utilizing the TV!) scheme for the Navier- 
Stokes equations is to difference the hyperbolic terms the same way as before, and then central 
difference the viscous terms. The final algorithm is the same as equation (3.1) except that 
the spatial central differencing of the viscous term is added to the right hand side of equation 
(.U . I he numerical solution shown below illustrates that this algorithm produces a fairly good 

solu ion for the case of a R A102822 airfoil calculation. A treatment for time-accurate calculations 
can be found in reference | 2 ], 

Numerical Results 

Generali y, for inviscid calculations, upwind TV I) schemes produce sharper shocks than sym¬ 
metric I VI) schemes |7|. For the current two methods ( 2 . 6 ) and (2.7), this seems to be not the 
case. I he two methods appeared to produce almost identical results for flow field conditions 
ranging from subcntical to transonic and supersonic. Figures 1 and 2 show the comparison of 
equation ( 2 . 6 ) with (2.7b) for two inviscid steady-state airfoil calculations. The advantages of 
symmetric TVI) schemes arc that they require less computational effort-and provide a more 
natural way of extending the scheme to two and three-dimensional problems. 

Figures 3 and 4 show an inviscid comparison of the symmetric TVI) scheme with the widely 
distributed computer code ARC 2 I), version 150 ( 9 ). The freestream Mach numbers are A/ c , | }> 
and 1.8, and the angle of attack is a - 7. The pressure coefficient distributions (not shown) are 
ldcntica between the two methods and yet the flow field appears very different. The symmetric 
1 VI) scheme gives a very well-ordered flow structure and can still capture the shocks with a 
coarse grid, especially near the trailing edge of the airfoil. On the other hand, the ARG2I) 
code did rather poorly. The same problem was studied for the upwind TVI) scheme and the 
results and convergence rates were found to be almost identical to those for the symmetric TVI) 

SC | om!™ SI<il,al ° f 10 12 Can be reachcd at arou,,d 400-506 steps. ARC21), however, required 
only 200-300 steps to converge to the same residual. The ARC2D, version 150 computer code is 

based on the Beam and Warming ADI algorithm | 10 | but uses a mixture of second ami fourth- 
order numerical dissipation terms. Figure 5 is an example of the viscous case for the RAI02822 
air o, using the upwind TVI) scheme. The thin layer Navier-Stokes equations with the algebraic, 
ur mlence model of Baldwin and Lomax | 12 | are used and the transition is fixed at 3 % chord 

I he overall agreement with experiments is quite good. The i 2 -norm residual of 10 ‘ 7 can be 
reached in around 900 steps. 

Concluding Remarks 

Both the symmetric and upwind TVI) schemes are designed to capture shock waves accurately 
while not exhibiting the spurious oscillation associated with the more classical second-order 
schemes. Numerical experiments with the TVI) schemes on problems containing no shock (ill 
show that there is no advantage of the TVI) scheme over the conventional Lax-Wcndroff type 
scheme. Numerical experiments with various inviscid airfoil calculations show that for weak to 
moderate shocks the main advantage of the TVD schemes is that one can capture the shocks 
m 1-2 grid points |5|. The flow field away from the shock looks very much like the classical 
second-order central diflerence methods. But as one increases the shock strength, especially in 
supersonic flow, 'I'VI) schemes provide superior flow field solutions even with a very coarse vet. 
non-clustering grid near shock waves. Numerical experiments also indicate that the symmetric 
and upwind i’VI) schemes are also applicable for viscous calculations. The current study further 
shows that the symmetric TVI) scheme is just as accurate as the upwind TVD scheme while 
requiring less computational effort than the upwind TVI) scheme. 
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